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The Objective of the Test: 

PURPOSE 
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Purpose 


• NASA still depends heavily on the attenuation methods of the Pyrotechnic Design 
Guidelines Manual for preliminary pyroshock environment estimation. 

• Project Goal: Understand and quantify how various composite panel properties impact the 
composite panel's response to a pyroshock environment. 



tSUiXtHOS TO? HONEYCOMB STRilCTUE - . 



Kacena, W. J., McGrath, M. B., & Rader, W. P. (1970). 
Pyrotechnic Shock Design Guidelines Manual. Denver: 
NAS5-15208. 
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Purpose Cont'd... % 


• Project Process: Use design of experiments (DOE) techniques to quantify differences in 
effects and variability in responses due to changes in input factors. 

• Challenge: 

• Shock environments can be very difficult to quantify. The most common methods are... 

• Time history based: shock pulses (half-sine pulse, etc.), wavelet reconstruction,... 

• Response spectrums: SRS, pseudovelocity (PV) 

• Analysis methods team starting with single value inputs; shock tends to be 
characterized with a spectrum. 

• Response spectrums for pyroshock tend to show similar trends (slope, frequency break 
point, plateau), but enveloping response spectrums can be very subjective. 
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Purpose Cont'd... £ 


•Goal 1: Eliminate human subjectivity by automating 

• Post-Processing of acceleration time history data 

• Use first 20 ms 

• Remove bit-error 

• Detrend 

• Enveloping of Shock Response Spectrum 

• Enveloping of Pseudovelocity 

• Goal 2: Determine and calculate single value inputs that can characterize the shock 
environment 

• SRS Envelope Parameters 

• PV Envelope Parameters 
•Temporal Moments 

• Spectral Moments 

• Goal 3: Analyze additional spectrums used for characterizing shock 

• Fourier spectrum 

• Energy Spectral Density 

• Time-Frequency spectrum 

• Goal 4: Utilize statistical processes to evaluate the data 

• Isolate non-essential parameters 

• Develop scaling method for composite structures 
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Collecting the needed data: 

TEST SETUP & DESIGN 
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Group I - 

Monolithic 

Composite 

Thickness 

(inch) 

Ply and Orientation 

Fill 

LSC" Core 
Load 
(gpf) 

Number 
of Tests 


Thin 

0 degree. 

unidirectional, tape 

N/A, solid 

10 

4 

Thin 

Quasi Isotropic, fabric 

N/A, solid 

10 

4 

Thin 

Quasi Isotropic, tape 

N/A, solid 

22 

4 

Thick 

Quasi Isotropic, tape 

N/A, solid 

10 

4 

Thick 

Quasi Isotropic, fabric 

N/A, solid 

22 

4 

Group II - 

Filled 

Composite 

Thickness 

Ply Orientation 

Fill 

LSC Core 
Load 



8 Ply Fill 

Quasi Isotropic, tape 

Al 

Honeycomb 

10 

4 

8 Ply Fill 

Quasi Isotropic, tape 

Rohacell 

Foam 

10 

4 
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Test Outline 



Group I - Solid Composite Panels 

Test 

Number 

Material 

Panel 

Thickness 

Ply 

Orientation 

Type 

LSC' 

Core 

Load 

1 

Composite, 

IM7/TC350 

Thin 

Fabric 

0-Deg. 18 ply 

Solid 

10 

2 

Composite, 
IM7/TC3 50 

Thin 

Fabric 

+45%45°, 0° (2x), 
+45%45°, 

90° (2x), 18 ply 

Solid 

10 

3 

Composite, 

IM7/TC350 

Thick 

Tape 

+45°/-45°, 0° (2x), 
+45°/-45°, 

90° (2x), 54 ply 

Solid 

10 

4 

Composite, 

IM7/TC350 

Thick 

Fabric 

+45°/-45°, 0° (2x), 
+45°/-45°, 

90° (2x), 27 ply 

Solid 

22 

5 

Composite, 

IM7/TC350 

Thin 

Tape 

+45%45°, 0° (2x), 
+45°/-45°, 

90° (2x), 38 ply 

Solid 

22 

Analyze results 

rom Tests 1-5 and re-plan as necessary 

6 

Composite, 

IM7/TC350 

Thin 

Fabric 

0-Deg. 18 ply 

Solid 

22 

7 

Composite, 

IM7/TC350 

Thin 

Fabric 

+45%45°, 0° (2x), 
+45°/-45°, 

90°, 18 ply 

Solid 

22 

8 

Composite, 

IM7/TC350 

Thick 

Fabric 

+45°/-45°, 0° (2x), 
+45°/-45°, 

90° (2x), 27 ply 

Solid 

10 

9 

Composite, 

IM7/TC350 

Thick 

Tape 

+45°/-45°, 0° (2x), 
+45°/-45°, 

90° (2x), 54 ply 

Solid 

22 

10 

Composite, 
IM7/TC3 50 

Thin 

Tape 

+45°/-45°, 0° (2x), 
+45°/-45°, 

90° (2x), 38 ply 

Solid 

10 
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Group II - Sandwich Composite Panels 

Panel 

Number 

Test 

Order 

Number 

Material 

Panel 

Thickness 

Fill Ply 

Orientation 

Type 

LSC’ Core 
Load 

11 

1 

Composite. 

IM7/TC350 

8 Ply 
Fill 

Al 

Honeycomb 
& Tape 

90°/+45 o /-45 o /0°/0°/ 

-45°/+45°/90°, 

8 ply both faces 

Sandwich 

10 

12 

4 

Composite. 

IM7/TC350 

8 Ply 
Fill 

Al 

Honeycomb 
& Tape 

90 0 /+45 0 /-45 0 /0°/0 0 / 

-45°/+45°/90°, 

8 ply both faces 

Sandwich 

22 

13 

3 

Composite. 

IM7/TC350 

8 Ply 
Fill 

Rohacell 
Foam &Tape 

907+45%45°/0°/0°/ 

-457+45790°, 

8 ply both faces 

Sandwich 

10 

14 

2 

Composite. 

IM7/TC350 

8 Ply 
Fill 

Rohacell 
Foam &Tape 

907+45%45°/0°/0°/ 

-45°/+45°/90°, 

8 ply both faces 

Sandwich 

22 

15 

7 

Composite. 

IM7/TC350 

8 Ply 
Fill 

Al 

Honeycomb 
& Tape 

90 o /+45 o /-4570 o /0°/ 

-457+45790°, 

8 ply both faces 

Sandwich 

10 

16 

5 

Composite, 

IM7/TC350 

8 Ply 
Fill 

Al 

Honeycomb 
& Tape 

907+457-4570707 

-457+45790°, 

8 ply both faces 

Sandwich 

22 

17 

8 

Composite, 

IM7/TC350 

8 Ply 
Fill 

Rohacell 
Foam &Tape 

907+457-457070°/ 

-457+45790°, 

8 ply both faces 

Sandwich 

10 

18 

6 

Composite, 

IM7/TC350 

8 Ply 
Fill 

Rohacell 
Foam &Tape 

907+457-4570707 

-457+45790°, 

8 ply both faces 

Sandwich 

22 


Evaluate results from tests 11-18 and determine test panel configuration for tests 19-28 
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Test Article 




Test Panel 
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Ceiling 



Floor 
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Characterizing the Test Data: 


SINGLE VALUE INPUTS 
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Single Value Inputs - Shock 

Response 



• A shock response spectrum can often 
be enveloped knowing three parameters 
(assuming a frequency range of 100 to 
10,000 Hz): 

• Frequency break point (Hz) 

• Max Peak Accel or plateau value (G) 

• Slope (dB/oct) 
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• A pseudovelocity spectrum can be 
enveloped knowing three 
parameters (assuming a frequency 
range of 100 to 10,000 Hz): 

• Mean Pseudovelocity (ips) 

• Max Peak Accel or plateau 

value (G) 

• Slope (dB/oct) 


Pseudovelocitv 


Mean PV 



Naturalf [Hz] 
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Single Value Inputs - Temporal 
& Spectral Moments 



Single Value Inputs: Temporal and Spectral Moments 

Temporal Moment Calculated Values 

Spectral Moment Calculated Values 

Temporal Energy 

Spectral Energy 

Temporal Mean 

One Sided Spectral Mean 

Temporal Variance 

One Sided Spectral Variance 

Root-Mean-Square Duration 

One Sided RMS Bandwidth 

Variance Normalization 

Variance Normalization 

Temporal Skewness 

One Sided Spectral Skewness 

Temporal Kurtosis 

One Sided Spectral Kurtosis 

Root Energy Amplitude 
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Reducing human error: 

AUTOMATING THE PROCESS 
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SRS Enveloping - General ^ 

IV ^ 

Algorithm ^ 


• Read in SRS plot points. 

• Determine a temporary frequency break 
point (Lowest frequency that is one standard 
deviation below the highest peak value). 

• Calculate the slope of the data points from 
100Hz to the temporary frequency break point 
(Least Squares Fit). 

• For the plateau, calculate the mean value of 
the data points from the temporary frequency 
break point to 10,000Flz. 

• Calculate where the sloped line and the 
plateau intersect; this is the new frequency 
break point. 

• Create the sloped portion of the curve from 
lOOFIz to the new frequency break point. 

• Create the plateau from the new frequency 
break point to 10,000Hz. 

• Output the frequency break point, the slope, 
and the plateau values. 


NEST1-01 ,048 - Ch 3 
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PV Enveloping - General ^ 

IV ^ 

Algorithm ^ 


• Read in PV plot points. 

• Determine break points (lowest and highest 
frequencies that correspond to peaks that are 
two standard deviations below the highest 
peak.) 

• Calculate the slope of the max displacement 
line. 

• Calculate the mean of the max 
pseudovelocity line. 

• Calculate the slope of the max acceleration 
line. 

• Calculate the two break frequencies by 
calculating where the max displacement and 
max acceleration line intersect the max 
pseudovelocity line. 

• Create the enveloping curves. 

• Output the max displacement value, the 
pseudovelocity, and the max acceleration 
value. 


NEST1-01 1 048-Ch3 
Pseudovelocity Envelope: Q=10 



Frequency, Hz 
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Post Processing Automation - 
General Algorithm 


NASA 


• Import time history data. 

• Identify the beginning of the shock pulse in the time history. 

• Remove the bit error from the time history. 

• Take the first 20msec of the shock pulse. 

• Detrend (linear) the shock pulse. 

• Perform calculations on post-processed data. 

• Time history plots 

• Shock Response Spectrum 

• Pseudovelocity Response Spectrum 

• Temporal Moments 

• Spectral Moments 

• Fourier Spectrum 

• Generalized Harmonic Wavelet Transform 

• Print and plot results for statistical evaluation. 
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Searching for Trends: 

PRELIMINARY STATISTICAL 
ANALYSIS 
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Max Acceleration Stats Analysis: 
Adequate Model Significant Factors 


NASA 


• Distance (top 
graph) has an 
overall effect of 
decreasing from 6 
to 42 inches, then 
leveling out 


Maximum Acceleration, Overall Means 


c 

o 

« 

k. 

_a> 

v 

o 

o 

< 


E 

3 

E 


x 

(0 



Distance, inches 


• The attenuation 
appears steeper 
and may bottom 
out at a lower 
value for Thicker 
panel than for 
Thinner (bottom 
graph) 


Maximum Acceleration, Panel Thickness 

—♦—Max Acceleration, 
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Max Acceleration Stats Analysis: 
Adequate Model Significant Factors 


NASA 


• LSC Core Load 
may or may not 
be a significant 
factor (top) 

• Ply does not 
appear to be a 
significant driver 
of Maximum 
Acceleration, 
given this model 
(bottom) 


Maximum Acceleration, LSC Core Load 



E 

D 

E 

x 

g 1 1 1 1 1 1 — 

“ 6 18 30 42 54 66 

Distance, inches 

Maximum Acceleration, Ply 


—♦—Max Acceleration, 
Ply = Fabric 
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Maximum Acceleration 



Max Acceleration Stats Analysis: 
Graph of Data with Full Adequate Model 



N'A 


Maximum Acceleration: Data and Overall Expected 

Values by Distance 


t=thin 

T=thick 


o 

o 

□ 


o 


X 

A 

Data, Panel 2 

-a ' t .Fabric, 10 


□ 

Data, Panel 3 

— T .Tape, 10 


- 

Data, Panel 4 

T .Fabric, 22 


+ 

Data, Panel 5 

t .Tape, 22 

o 

X 

Data, Panel 7 

• °x>* t .Fabric, 22 


o 

Data, Panel 8 

— o— T .Fabric, 10 


O x 

Data, Panel 9 

—x— T .Tape, 22 


0 

Data, Panel 10 

-o- t .Tape, 10 



n 1 1 1 1 r~ 

6 18 30 42 54 66 

Distance from Panel Edge, inches 
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Temporal Energy: Visual Analysis to>a 






In general, there 
appears to be a 
decrease in Temp 
Energy with Distance 

Tape may show a 
larger response than 
Fabric - or we may be 
seeing noise 

Top/ Bottom seem to 
replicate each other 
well 

0 and 45/90 appear to 
replicate one another 


07/19/2013 4:02 PM 

Data Table=130705 131 1 SBU Composite ShockProcessed 
Results Tests 1 - 12 

Graph Builder 


C 

LU 


Temp Energy vs. Ply 
Panel Thickness 


Thin 


Thick 



c/> 

o 

o 

o 


o 

D> 

CL 


Where(24 rows excluded) 


Distance 
— 1 
— 2 

3 

4 

5 

6 
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• Thicker panels appear to 
have lower Temp Energy 
than Thin 

— The effect is very noisy 
- may be difficult to 
pick out in quantitative 
analysis OR may not 
exist 

- Panel 2 (Thinner Fabric 
10 gpf) might have 
given low values - 
factor combination not 
on retest list, though 
Panel 2 is 

• Panel 5 (Thin Tape 22 
gpf) seems to be noisy 
close-in 

- Likely not of great 
concern to analysis, 
but could be useful to 
know why 


07/19/2013 4:02 PM 

Data Table=130705 131 1 SBU Composite ShockProcessed 
Results Tests 1 - 12 

Graph Builder 


Temp Energy vs. Panel Thickness 
LSC Core Load 

10 22 



Distance 
— 1 
— 2 

3 

4 

5 

6 
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• Visual analysis: 

- Usually, Core Load 
doesn't make a 
difference 

- Usually, Tape gives 
a higher response 
than Fabric 

- Usually, Panel 
Thickness doesn't 
matter 

- BUT Thin Fabric is 
quite sensitive to 
Core Load (3-way 
interaction) 

• BUT data is equivocal 


07/09/2013 1:45 PM 

Data Table=1307081346 SBU Composite ShockProcessed 
Results Tests 1 - 10 

Graph Builder 


£ 


Mean PV vs. LSC Core Load 
Ply 


Fabric Tape 


2 

2 

Is this apparent 
response due to 
when these tests 
were run? 

* 

I 1 

A 

• 

; • 

» 1 

+ 

+ 

+ 

+ 

+ 

1 t 


t 1 1 1 1 1 1 — — i 1 1 1 1 1 r 


■o 

s> 

9 

S 


Test 
• 1 
• 2 

• 3 

• 4 

• 5 

• 6 

• 7 

• 8 

• 9 

• 10 


LSC Core Load 
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• Mean PV by Test 
Number 

- Range of 
variability - 
suggests 
measurement is 
noisy in particular 
instances 

— What are the 
instances? 

• In particular, Test 9 
has several wild 
points 


07/09/2013 1:45 PM 

Data Table=130708 1346 SBU Composite Shock Processed 
R esults Tests 1-10 

Graph Builder 


Mean PV vs. Test 
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• Max Accel generally decreased with increasing Distance 

— Some panels showed a knee at ~18", with maximum values 
there for a few panels 

• Little difference seen between Fabric and Tape panels 

• Thin panels gave higher Max Accel responses than Thick 
Panels over a longer distance 

— Thick panels' Max Accel tended to fall off steeply before 18" 

• 22 gpf Core Load tended to result in lower Max Accel 
values at 6" than 10 gpf 

— Thin panels' Max Accel values increased to a maximum at 18" 
Distance 
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• Top accelerometer TE values are considerably and 
fairly repeatably higher than Bottom 

• TE generally decreases with decreasing Distance 

— Some panels' curves show a knee at ~30" 
Distance, with a few even showing a maximum 
there 


• Thinner Thickness results in higher TE than Thicker 

• Variability increases with increasing TE 

— Analyzed logarithm of TE - analysis was of little 
use without this transformation 
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Mean PV: Bottom Lines 



• There are a number of wild points that seriously inhibit 
quantitative analysis. 

— Recommend looking at these points in the data to see if 
there is something driving this. 

— Test 9 exhibits particularly high variability. 

• In visual analysis, after attempting to identify and 
disregard these wild points, the following conclusions 
might be made: 

— Core load and Panel Thickness usually don't make much 
difference. 

— Fabric usually gives lower Mean PV than Tape. 

— BUT Fabric at Thinner Thickness and 22 gr Core Load 
appears to have higher Mean PV. 

- 0-degree panels may respond similarly to 45/90 panels. 
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• Evaluate noise in the data. 

• Check if algorithms are too sensitive. 

• Review time histories and any other factors that might 
explain noise and outliers. 

• Search for trends in data. 

• Understand the physical meaning of parameters. 

• Some pseudovelocity results are counter to expectations. 

• Consider retesting panels after non-destructive evaluation of 
panels. 

• Complete all test series. 

• Develop useful tool for a composite panel's response to a shock 
environment. 

• Identify non-significant parameters of composite panels. 

• Develop scaling methods for composite panels, if possible. 
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Questions? 
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Back up 
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Single Value Inputs - Temporal 

Moments 


Accel = Accelerati onTimeHist oryi^A) 

•Temporal Energy ,Tnrg = ^ Accel 2 * At 


Normalization, Anorm = 


Accel 

Tnrg 


Temporal Mean, Tmean =^t* Anorm * At 


•Standardization, tO = t -Tmean 


• Temporal Variance, T var = ^ 

• RMS Duration, D = Vt var 

• Variance Normalization, tO = 


tO 2 * Accel 2 * At 
tO 

Z 3 ^ 

1 0 *Accel^*At 


• Root Energy Amplitude, 7m? 



• Temporal Kurtosis, Tkurt = * Accef * At 


Hacker, J. (2012, May 30). Index of ula time-frequency motlab scripts. 
Retrieved from ftp://shockwg@drop. aero.org/ 
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Single Value Inputs - Spectral 

Moments 


Accel = Accelerati onTimeHist ory{\) 
x (/) = fftshift (fft(Accel)) * At 
XX(f) = \X(f)\ 2 

• Spectral Energy, Fnrg = ^ XX (/) * A / 

• Normalization ,Xnorm = ^ill 

Fnrg 


• One Sided Spectral Mean, Fmean 1 - ^ (j/| * XX (/))* A/ 

• Two Sided Spectral Mean, Fmean 2 = £(/*AX(/))*A/ 

• Standardization, /l = / _ Fmean 1 

/2 = / 

• One Sided Spectral Variance, F var 1 = £(/l 2 *AX(/))*A/' 

• Two Sided Spectral Variance,^ 7 var 2 = £(/2 2 *AT(/))*A/ 

Hacker, J. (2012, May 30). Index of ula time-frequency matlob scripts. 
Retrieved from ftp://shockwg@drop.aero.org7 
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Single Value Inputs - Spectral 
Moments Cont'd... 


• One Sided RMS Bandwidth, B\ = y/F varl 

• Two Sided RMS Bandwidth, B2 = V^var2 

• Variance Normalization, = f 2 = 12l 

J B\ B\ 


• One Sided Spectral Skewness, Fskewl = ^(|/l| 3 *AX(/))* A/ 

• Two Sided Spectral Skewness, Fskewl = X(/2 3 *AX(/))*A/ 

• One Sided Spectral Kurtosis, Fkurtl = X(/1 4 *AX(/))*A/ 

• Two Sided Spectral Kurtosis, Fkurtl = X(/2 4 *XY(/))*A/ 



Hacker, J. (2012, May 30). Index of ula time-frequency motlob scripts. 
Retrieved from ftp://shockwg@drop.aero.org7 
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Results - Output Figures 





0 0 005 0.01 0015 0 02 0 025 

Time (s) 



0 0 005 0.01 0.015 0 02 0.025 

Time (») 


| 

t 

I 


NESTI-02,048 - Ch 1 Fourier Spectrum 



Frequency [Hz] 


NEStT^oiaTon 

SRS Envelope Companion Q=10 



NEST1-02,04S-Ch1 
Pseudoveloclty Envelope; Q=10 
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Results - Output Figures 

Cont'd.... 




Time (s) 



x 10 


Temporal Density - Amplitude ((in/s 2 ) 2 ) 



Amplitude ((in/s 2 )J:-$£Rz) 


Time (s) 


Hacker, J. (2012, May 30). Index of ula time-frequency motlob scripts. Retrieved 
from ftp://shockwg@drop. aero.org/ 
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SRS Enveloping - Code ~W 


Contents 

■ Determine the intermediate frequency break point of the slope and plateau 
« Find envelope line properties 

• Create envelope curve 

■ Plot results 

function [ env ] = srs_envelope2 ( srs, varargin ) 


%srs_envelope Envelopes an SRS. Creates a plot and returns the values of 
Ithe break point, slope, and plateau value 
% SRS i3 a two-column matrix of frequency and G values 

f=srs (: , 1) ; 
g=srs (: ,2) ; 


Error using srs_envelope2 (line 6) 
Not enough input arguments. 


Determine the intermediate frequency break point of the slope and plateau 

[pics, locs]=findpeaks (srs ( : , 2) ) ; 
s=std(plcs) ; 

pl=zeros (length (pks) ,1) ; 
ii=i; 

for ii=l : length (pks ) 

if pk3(ii)>=(max(pks)-l*3) 
pl(jj)=loc3(ii) ; 

33=33+1; 

end 

end 

% remove zeros from pi 
pi (pl=0) = [] ; 

I now use the first point as the initial frequency break point 
fbi=pl (1) ; 
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SRS Enveloping - Code 


N 



Find envelope line properties 


find line properties from OHzto the initial breakpoint 


mdl=LinearModel. fit (loglO (f(l:fbi) ) , loglO (g(l:fbi) ) ) ; 

% y=b*x A N 

N=mdl. Coefficients (2, 1} ; lintercept 
b=10 A mdl.Coefficients(l,l}; Islope for exponent 

% Find mean value for the plateau 

% Note: mean was chosen rather than max because it will take into account 
I all of the values in the plateau range. 
yp=mean (g ( fbi : end) ) ; 

% Calculate break point 
fb=(yp/b) A (l/N); 


Create envelope curve 

Determine counter for envelopes 


ss=l; 

while f(ss)<fb 
s3=ss+l; 
if ss>length(f) 
break 

end 

end 

ss=ss-l; 

I Sloped part of line 
fslope=zeros (ss, 1) ; 
yslope=zeros (ss, 1) ; 
for gg=l:ss; 

fslope {gg)=f (gg) ; 
yslope (gg) =b*f (gg) A N; 

end 

1 Remove excess zeros 
fslope (fslope=0) = [] ; 
yslope (yslope ==0 ) = [] ; 

% Calculate the slope in dB/Oct 

dbOct=2*loglO (20) *logl0 (yslope (end) /yslope (1) ) /loglO (fslope (end) /fslope (1) ) ; 

I Plateau 

cc=l; 

fplateau=zeros (length (f) -ss, 1) ; 
yplateau=zeros (length (f ) -33, 1) ; 
for hh=ss: length (f) 

fplateau (cc) =f (hh) ; 
yp 1 a t e au ( cc ) =yp ; 
cc=cc+l; 

end 

% Remove excess zeros 
fplateau (fplateau=0 ) = [ ] ; 
yplateau (yplateau=0 ) = [ ] ; 
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Plot results 

loglog(f,g, 'LineWidth' ,2) 
hold on 

loglog(fslope,yslope, 'r' f 'LineWidth' , 2) 
loglog(fplateau,yplateau, 'g', 'LineWidth' ,2) 
grid on 

title ( [varargin, ' SRS Envelope Comparison; Q=10 ' ] ) ; 

legend (' Shock Response Spectrum' , 'Constant Velocity Line' ,' Plateau' , 'Location' , 'Best' ) ; 
xlafcel ( 'Natural Frequency (Hz)') 
ylafcelC Peak Accel (g) ' ) 
hold off 

I Output envelope properties 

env=[fb,dbOct,yp] ; 


end 
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SRS Enveloping Algorithm Limitations 


N 


1 
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SRS Enveloping Algorithm Limitations 


N 


1 


Stage 1 Ignition Event 



Stage 1 Ignition Event 
SRS Envelope Comparison: Q=10 



Natural Frequency (Hz) 
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PV Enveloping - Code 


Contents 

■ Determine the frequency break points 

■ Create curve envelope 

■ Create figure 

■ Calculate PV parameters 



function [ env ] = pv_envelope2 ( pseudo, varargin ) 


lpv_envelope2 Envelopes a pseudovelocity plot and outputs a figure of the 
%enveloped curve and the primary properties of the curve. 


Determine the frequency break points 

[pks, Iocs] =findpeaks (pseudo ( : , 2) ) ; 
s=std(pks) ; 

pl=zeros (length (pks) , 1) ; 

33-1; 

for ii=l: length (pics) 

if pks (ii)>= (max (pics) -2*s) 
pl(jj)=locs(ii); 

33-33+1; 

end 

end 

% Remove zeros from pi 

pi (pl— o ) = [ ] ; 

I Now use the first and last points for your bounds 
k=[l,pl(l) ,pl(end) , length (pseudo (:, 1) ) ] ; 

I Find max displacement line properties 

mdl_l=LinearMcdel.fit(loglO(p3eudo(k(l) :k(2) ,1) ) , loglO (pseudo (k(l) :k(2) , 2 ))); 

% y=b*x A N 

N_l=mdl_l . Coefficients {2,1}; lintercept 
b_l=10 A mdl_l.Coefficients{l,l}; Islope for exponent 

1 Find mean value max pseudovelocity plateau 

I Note: mean was chosen rather than max because it will take into account 
% all of the values in the plateau range. 
yp=mean (pseudo (k (2) :k(3) ,2) ) ; 

% Calculate first break point 
fb_l= (yp/b_l) A (1/N_1) ; 

% Find max acceleration line properties 

mdl_2=LinearModel. fit (loglO (pseudo (k (3) :k(4) ,1) ), loglO (pseudo (k (3) :k(4) ,2) ) ) ; 
I y=b*x A N 

N_2=mdl_2. Coefficients {2, 1} ; % intercept 
b_2=10 A mdl_2.Coefficients{l, 1} ; Islope for exponent 

% Calculate second break point 
fb_2= (yp/b_2) A (1/N_2) ; 
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PV Enveloping - Code 


Create curve envelope 

Determine counters for envelopes 



ss_l=l; 

while pseudo (ss_l, l)<fh_l 

33_1=3S_1+1; 

if ss_l>length (pseudo (: , 1) ) 
break 

end 

end 

33_1=3S_1-1; 

ss_2=l; 

while pseudo (S3_2, l)<fb_2 
33_2=ss_2+1; 

if ss_2>length (pseudo (:, 1) ) 
break 

end 

end 

ss_2=ss_2-l; 

% Sloped part of max displacement line 

mdispl=zeros (ss_l, 2) ; 
for gg_l=l:ss_l; 

mdispl (gg_l, l)=pseudo (gg_l, 1) ;%frequency 

mdispl (gg_l, 2)=b_l*pseudo (gg_l, 1) ~N_1; ipseudovelocity 

end 

I Pseudovelocity plateau 
cc=l; 

pvplateau=zeros (ss_2-ss_2, 2) ; 
for hh=ss_l:s3_2 

pvplateau(cc, l)=pseudo (hh, 1) ; %freguency 
pvplateau (cc, 2) =yp; %p3eudovelocity 
cc=cc+l; 

end 

% Sloped part of max acceleration line 

maccel=zeros (length (pseudo ( : , 1) ) -ss_2, 2) ; 
for gg_2=ss_2 : length (pseudo ( : , 1) ) ; 

maccel (gg_2, l)=pseudo (gg_2, 1) ; %frequency 

maccel (gg_2, 2) =b_2 ‘pseudo (gg_2, 1) ^N_2; %pseudovelocity 

end 
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PV Enveloping - Code "W 


Create figure 

loglog (pseudo ( : , 1) , pseudo (: ,2) , 'b* f 'LineWidth' ,2) 
hold on 

loglog (mdispl (: f 1) , mdispl ( : , 2) , 'm' , 'LineWidth' ,2) 
loglog (pvplateau ( : , 1) ,pvplateau( : , 2 ) , ' g' , 'LineWidth' , 2) 
loglog (raaccel(: , 1) ,maccel (: , 2 ) , ' r' , 'LineWidth' ,2) 
loglog (pseudo (le, 1) , pseudo (k, 2) , ' k* ' , 'LineWidth' ,2) 
hold on 
FourcpDP 
hold off 

title ( [varargin, ' Pseudovelocity Envelope; Q=10']) 

legend (' Fseudovelocity' , 'Max Displacement Line', 'Mean PV','Max Accel Line 'Border Points 'Location' , 'Best' ) ; 


Calculate PV parameters 

Max value of max displacement line 


maxdisp=max (mdispl (: , 2) ./ (2*pi () *mdispl (: , 1) ) ) ; 

I Max value of max acceleration line 

maxg=max (maccel ( : , 2 ) . *maccel ( : , 1) *2*pi() /386. 1) ; 

env= [maxdisp, yp,maxg] ; 


end 


2013 - 11-07 


NASA: MSFC 


48 






Post Processing Automation 


% Choose a name that the data and outputs will be saved under. 
test_set=' Group- lJTest-2' ; 


Prepare Accel Time Data for Post Processing or load existing processed data 

First, checkto see if the data has already been saved in a mat file of the name test_set. 



check=dir ( [test_set, * .mat' ] ) ; 
if isempty (check) =1 
check=( [] ) ; 
check. name=* false' ; 

end 

if strcrapi (check. name, [test_set f ’ .mat ' ] ) =1 
load ( check . name ) 

else % If data has not already been saved, import the data 
% Import .csv file information 
info=dir ( ' * . csv' ) ; 
test=struct ( [] ) ; 
for ff=l: length (info) 

test (ff) .name=regexprep (info (ff, 1) .name, ' - Time Data. csv' ; 
test(ff) .data=csvread(info (ff , 1) .name, 5,0); 

end 

clear info ff 

% Find where shock pulse begins by locating the first data point that is 
% Greater than the bit error (assume first 10,000 data points). 

counter=ones ( 1 , length ( tes t ) ) ; 
if length (test (1, 1) .data (: , 1) ) >=100000; 
for ff=l: length (test) 

for jj=l: length ( test (l,ff) .data(:,l) ) 

biterror=max(abs (test (1, ff ) .data (1:10000, 2) ) ) ; 
if test(l,ff) .data(jj,2)<=biterror 
counter (ff) =counter (f f ) +1; 

else 

break 

end 

end 

if counter (ff) >=100000 
counter (ff) =60000; 

end 

end 

% Remove bit error from data then take detrended 20ms after pulse 
for ff=l: length (test) 

test (1, ff ) .data=[test (1, ff ) .data ( : , 1) , test (1, ff ) .data(: ,2) -mean (test (1, ff ) . data (1: counter (ff) , 2) ) ] ; 

test(l,ff) .data=[test (1, ff) .data (counter (ff) : counter (ff) +20000, 1) ,detrend(test (1, ff ) .data (counter (ff) : counter (ff) +20000, 2) ) ] ; 
% Shift time so that it starts at zero 
dt=test (l,ff) . data (2,1) -test (l,ff) .data (1,1) ; 

N=length(test(l,ff) .data(:,l) ) ; 

test(l,ff) .data=[ (0 :N— 1) * *dt, test (1, ff ) .data(:,2)]; 
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/"* Post Processing Automation Cont'd.. 




lata is now ready to be post processed. 

Save figure of time histories (accel, veloc, disp) 


for ff=l: length (test) 

timehistplot (test (1, ff) .data, test (1, ff) .name) ; 

set (gcf, 'Units 'normalized' , 'OuterPosition' , [0 0 1 1]); 

saveas (gcf, [test (1, ff ) .name, '_time' ] , 'jpeg’ ) ; 

close; 

end 

clear ff 
elf 

% Calculate SRS and SRS envelop and save figures 

for ff=l: length (test) 

test(ff) .srs=srsfunc(test(l,ff) .data) ; 

test (ff ) . srs_env_prop=srs_envelope2 (test (1, ff ) . srs, test (1, ff ) .name) ; 
set(gcf, 'Units', 'normalized', 'OuterPosition', [0 0 1 1]); 
saveas (gcf, [test (1, ff ) .name, '_srs_envelop' ] , ' jpeg' ) ; 
close 

end 

clear ff 

elf 

close 

% Calculate PV 

for ff=l: length (test) 

test(ff) .pv=pvssmax(test(l,ff) .data) ; 

test (ff ) .pv_env_prop=pv_envelope2 (test (1, ff) .pv, test (1, ff) .name) ; 
set(gcf, 'Units' , 'normalized', 'OuterPosition', [0 0 1 1]); 
saveas (gcf, [test(l, ff) .name, '_PV' ] , 'jpeg’ ) ; 
close 

end 

elf 

close 
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/"* Post Processing Automation Cont'd..,^ 


% Calculate fourier spectrum 

for ff=l: length (test) 

test(ff) . fr=fouriergab(test(l,ff) .data) ; 

semilogy (test (1, ff ) . fr (: , 1) , test (1, ff ) . fr (:,2) , 'LineWidth' , 2) ; 
grid on 

title ( [test (1, ff) .name, ' Fourier Spectrum']); 
xlabel ( ' Frequency [Hz ] ' ) 
ylabel ( 'Accel [in/sec y '2] ' ) 

set (gcf, 'Units' , 'normalized' , 'OuterPosition' , [0 0 1 1]); 
saveas (gcf, [test (1, ff) .name, '_Fourier' ] , 'jpeg' ) ; 
close 

end 

elf 

close 

% Calculate Temporal Moments 

for ff=l: length (test) 

test (ff ) . tm=tmoment_DP (test (1, ff ) .data ( : , 1) , test (1, ff ) .data ( : , 2) ) ; 

end 

% Calculate Spectral Moments 

for ff=l: length (test) 

test(ff) .sm=smoment_DP(test(l,ff) .data(:,l) ,test(l,ff) .data(:,2) ) ; 

end 

% Calculate Generalized Harmonic wavelet transform 
for ff=l: length (test) 

[test(ff) .S,test(ff) .F,test(ff) .H]=ghwt_DP (test (1, ff ) .data, test (1, ff) .name) ; 
set (gcf , 'Units' , 'normalized' , 'OuterPosition' , [0 0 1 1]); 
saveas (gcf, [test (1, ff) .name, '_ghwt ' ] , 'jpeg' ) ; 
close 

end 

elf 

close 
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Write parameters to a text file. 

fid=fopen( [test_set, '_Parameters.txt ' ] , 'w' ) ; 
fprintf (fid, ' ls\n' , test_set) ; 

fprintf (fid, ' Shock Response Spectrum Parameters'^' ) ; 

fprintf (fid, 'Output Name\t\t\t\t\t\t\t\t\t\t\tFrequency Break Point- (Hz) \t\tSlope- (dB/Oct) \t\tPeak- (g) \n' ) ; 
for ff=l: length (test) 

fprintf (fid, ' %s\t\t\t\t\t\t\t\t\t\t\t%6 . 2f\t\t\t\t\t\t\t\t\t%6 . 2f\t\t\t\t\t\t%6 . 2f \n' , test (1, ff ) . name, test (1, ff ) . srs_env_prop) ; 

end 

fprintf (fid, ' \nPseudovelocity Spectrum Parameters\n' ) ; 

fprintf (fid, 'Output Name\t\t\t\t\t\t\t\t\t\t\tMax Displacement- (in) \t\tMean Pseudovelocity- (ips) \t\tMax Accel- (g) \n' ) ; 
for ff=l: length (test) 

fprintf (fid, ' %s\t\t\t\t\t\t\t\t\t\t\t%5 . 3f\t\t\t\t\t\t\t\t\t%6 . 2f\t\t\t\t\t\t\t\ti6 .2f \n' , test (1, f f ) . name, test (1, f f ) .pv_env_prop) ; 

end 

fprintf (fid, ' \nTemporal Moments\n' ) ; 

fprintf (fid, 'Output Narae\t\t\t\t\t\t\t\t\t\t\tTemporal Energy\t\t\tTemporal Hean\t\t\tTemporal Variance\t\t\tTemporal Skewness\t\t\tTemporal Kurtosis\t\t\tRMS Duration\t\t\tRoot Energy Amplitude\n' ) ; 
for ff=l: length (test) 

fprintf (fid, ' Is\t\t\t\t\t\t\t%2.5e\t\t\t\t%2.5e\t\t\t\t\t%2.5e\t\t\t\t\t\t\tl2.5f\t\t\t\t\t\t\t\t%2.5f\t\t\t\t\t\t\t%2.5e\t\t\t\t\t%2.5e\n' , test (l,ff ) .name, test (1, ff) .tra(l, : ) ) ; 
end 

fprintf (fid, ' \nSpectral Moments\n' ) ; 

fprintf (fid, 'Output Narae\t\t\t\t\t\t\t\t\t\t\tSpectral Energy\t\t\tOne Sided Spectral Mean\t\t\tTwo Sided Spectral Mean\t\t\tOne Sided Spectral Variance\t\t\tTwo Sided Spectral Variance\t\t\tOne Sided Sp 
for ff=l: length (test) 

fprintf (fid, ' %s\t\t\t\t\t\t\t%2 . 5e\t\t\t\t\t\t%2 . 5e\t\t\t\t\t\t\t\t%3 . 4f \t\t\t\t\t\t\t\t\t%2 . 5e\t\t\t\t\t\t\t\t\t\t\t%2 . 5e\t\t\t\t\t\t\t\t\t\t\t%3 . 5f \t\t\t\t\t\t\t\t\t\t\t\t\t%3 . 5f \t\t\t\t\t\t\t\t\t\t\t\ 
end 

fclose (fid) ; 

save (test_set, 'test'); 
clearvars -except test 
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• Procedure 

— Look for potential model(s) using ANOVA on 
columnized dataset. 

— Examine model using repeated measures analysis. 

• Obtain model parameters. 

— Reanalyze using ANOVA on columnized dataset 
with model parameters from repeated measures 
analysis. 
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Max Accel: Visual Analysis ~W 


• Graph of 45/90 data 

- By input factor 
combination 

- Distance indicated by 
color 

- Log-scaled Max Accel 

• There appears to be 
decrease in Max Accel 
with increase in Panel 
Thickness 

- Decrease may be 
more pronounced 
with Tape than Fabric, 
but Tape may have a 
larger Max Accel at 
0.2 Panel Thickness 
than Fabric 

• There appears to be 
greater variability: 

- Closer-in Distances 

- 10 gpf Core load 


07/18/2013 1:09 PM 

DataTable=130717 1527 1 -12 repeat measures 

Graph Builder 

Max Accel 06 & 5 more vs. Panel Thickness 
Ply 

Fabric Tape 


E 

V) 

°a 



Max Accel 06 
Max Accel 18 
Max Accel 30 
Max Accel 42 
Max Accel 54 
Max Accel 66 


( n 
o 
o 


Where(4 rows excluded) 
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Max Accel: Visual Analysis ~W 


• Similar to previous 
graph, but showing 
individual points 

• Variability in Tape 22 
gpf now seems 
especially small 
compared to other 
factor settings 

• 0 degree appears to 
track 45/90 degree 
data 


07/15/2013 4:22 PM 

Data Table=130715 1607 Tests 1 - 12 

Graph Builder 


Max Accel vs. Panel Thickness 
Ply 

Fabric Tape 



Distance 
• 6 
• 18 

• 30 

• 42 

• 54 

• 66 

6 

18 

30 

42 

54 

66 
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Stats Analysis, Temporal Energy: 
Unconstrained Model 








Model with main effects and a 
number of interactions is 
adequate fit to data 

- Log TE modeled: says 
variability increases with 
increasing TE 

Significant, persistent 
difference between Top (darker 
lines) and Bottom (lighter) 

In general, Thinner thick panels 
give higher TE response than 
Thicker 

Panel 5 (green crosses, 0.2, 
Tape, 22) shows odd low values 
at close distance(s) 

- Max at 30" also seen in Panel 
9 (violet X's), possibly 2 (red 
triangles) and 7 (blue 
asterisks) 


>. 

U) 


c 

LU 


O 

Cl 

E 


Temporal Energy: Data and Modeled Expected Values by Distance 

Fabric, 10 B 
Tape, 10 B 
Fabric.22 B 
Tape, 22 B 
Fabric, 22 B 
Fabric, 10 B 
Tape, 22 B 
Tape, 10 B 



A Data, Panel 2 t Fabric.lOT t 

□ Data, Panel 3 m t Tape.lOT t 

- Data. Panel 4 t Fabric, 22 T t 

♦ Data, Panel 5 — t Tape, 22 T t 

* Data, Panel 7 t Fabric, 22 T t 

O Data, Panel 8 ■■ t Fabric.lOT t 

x Data, Panel 9 — x — t Tape, 22 T t 

o Data, Panel 10 —-<>-• t Tape.lOT t 
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Stats Analysis, Temporal Energy: 

Reduced Model 



• Constraining the 
model results in some 
ill fits, but on the 
whole it is probably 
still useful 


- Knee at 30" may not | 

represent some | 

panels' data - see | 

Panels 4 (yellow * 

dashes) and 8 
(indigo diamonds) 

• Largest responses 
now clearly result 
from Thinner 
Thickness 


Temporal Energy: Data and Modeled Expected Values by Distance 


X 


A Data. Panel 2 « — » — t Fabnc.lOT 

a Data. Panel 3 ' m t Tape.lOT 

- Data. Panel 4 - . r Fabnc,22T 

♦ Data. Panel 5 t Tape.22T 

x Data. Panel 7 • • « t Fabnc.22 T 

« Data. Panel 8 - ■«- t Fabnc.lOT 

x Data. Panel 9 — x— t Tape, 22 T 

o Data. Panell 0 t Tape.lOT 


x 


t Fabnc.lOB 
t Tape.lOB 
t Fabric,22 B 
t Tape.22B 
t Fabnc.22 B 
t Fabnc.lOB 
t Tape.22B 
t Tape.lOB 



18 30 42 

Distance from Panel Edge, inches 


54 


t=thin 

T=thick 
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Graphs of Mean PV Data NASA 


• Colored by 
Distance 

- Distance 
appears to 
have little 
effect in this 
measure 


07/09/2013 1:45 PM 

Data Table=l30708 1346 SBU Composite ShockProcessed 
Results Tests 1-10 

Graph Builder 

Mean PV vs. Test Distance 
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